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ABSTRACT 

We present a method, based on the correlation function of excursion sets above 
a given threshold, to test the Gaussianity of the CMB temperature fluctuations in 
the sky. In particular, this method can be applied to discriminate between standard 
inflationary scenarios and those producing non-Gaussianity such as topological defects. 
We have obtained the normalized correlation of excursion sets, including different levels 
of noise, for 2-point probability density functions constructed from the Gaussian, Xn 
and Laplace 1-point probability density functions in two different ways. Considering 
subdegree angular scales, we find that this method can distinguish between different 
distributions even if the corresponding marginal probability density functions and/or 
the radiation power spectra are the same. 
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1 INTRODUCTION 

The temperature field T{6,(j)), associated with the cosmic 
microwave background (CMB), is usually assumed to be a 
homogeneous and isotropic Gaussian random field, i.e. the 
n-point probability density function is a multivariate Gaus- 
sian characterized by the 2-point correlation function C{9) 
or equivalently the radiation power spectrum Ce. The form 
of the temperature field at recombination is related to the 
initial density field coming from quantum fluctuations dur- 
ing an infiationary era in the very early stages of the history 
of the Universe. However, other sources of density fluctu- 
ations (such as topological defects) may emerge at these 
early times and generate a non-Gaussian random fleld. In 
any case future high precision maps should be searched for 
any traces of non-Gaussianity since both foregrounds and 
systematic errors can leave non-Gaussian imprints as well. 

Whereas very accurate temperature maps can be ob- 
tained for the standard inflationary model (at the level of 
^ 1% in the radiation power spectrum), this is not the case 
for scenarios based on topological defects. In particular, the 
best known of these models is the one generated by cosmic 
strings (Turok 1996), which results in temperature maps ap- 
proximately Gaussian on scales that span from relative small 
(a few arcminutes according to Gott III et al. 1990) to large 



angular scales. Monopoles and textures generate maps which 
are also Gaussian at large angular scales but not at small 
scales (Shellard 1996). The texture and monopole maps show 
an asymmetric distribution with a positive tail decreasing 
more slowly than in the Gaussian case (Turok 1996). 

Taking into account the present uncertainty in char- 
acterizing the temperature distributions for some non- 
standard scenarios, in this paper we will explore several 
non- Gaussian models. We consider a Xn model for the tem- 
perature distribution function that satisfies the previous 
two properties, asymmetry and an exponential tail, and a 
Laplace model that represents a symmetric distribution with 
exponential tails. Moreover, we will also study a toy model 
which has the property of having a Gaussian marginal dis- 
tribution but a non-Gaussian 2-point probability density 
function, a property that appears in the cosmic string maps 
(Gott III et al. 1990). As it is shown in figures 12a,b of that 
paper, the temperature 1-point probability density function 
of cosmic strings is Gaussian on scales above a few arcmin- 
utes whereas the one of the temperature gradient shows 
a clear non- Gaussian behaviour. On the other hand, the 
patchy behaviour in a CMB map associated to such net- 
work can be represented by the mentioned toy model (see 
§2.2). 
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In the analysis of the 2D temperature maps, we will 
consider the correlation of the excursion sets (or regions) 
above a certain threshold. For the 3D matter density field, 
such correlations were introduced by Kaiser (1984) assum- 
ing a Gaussian 2-point probability density function to study 
the correlation function of rich clusters. He found an ampli- 
fication with respect to the correlation function associated 
with the underlying matter fluctuations. This was applied 
(Martinez-Gonzalez &: Sanz 1988) to test models of biased 
structure formation in several scenarios. Coles (1989) mod- 
eled the matter density field with a ID Gaussian, log-normal 
and a Xi distributions (the 2 non-Gaussian distributions 
were used in an attempt to allow for non-linear evolution) to 
study the correlation function of both peaks and excursion 
sets (in this last case, the Politzer & Wise 1984 approxima- 
tion was assumed). The excursion sets are easier to identify 
in the maps than the peaks and they carry very similar infor- 
mation at high thresholds. In fact, the number of maxima 
and excursion sets above a threshold coincides asymptoti- 
cally. 

For the 2D temperature field, non-Gaussian statistics 
have been applied by Coles & Barrow (1987) to obtain the 
mean size and frecuency of occurrence of the excursion sets 
above a given level for different random fields. Such prop- 
erties depend only on the 1-point probability density func- 
tion. On the other hand, Pompilio et al. (1995) used a mul- 
tifractal analysis of the temperature scans to explore the 
possibility of distinguishing between Gaussian fluctuations 
and non-Gaussian ones produced by a network of cosmic 
strings. Finally, Kogut et al. (1995) used the peak correla- 
tion function in the COBE-DMR 2-year sky maps to test 
for a class of non-Gaussian models. In the last paper, the 
temperature maps are generated using a spherical harmonic 
decomposition whose multipole components aim are drawn 
from Gaussian, log-normal and Xn distributions. They find 
that the 2-point correlation of peaks is a better discrimina- 
tor between Gaussian and non-Gaussian temperature fields 
than the genus and the 3-point correlation function. Kogut 
et al. (1996) also carried out a similar analysis with the DMR 
4- year sky maps. 

In a previous paper, we studied some geometrical prop- 
erties of the CMB maxima above a given level (mean num- 
ber, curvature and ellipticity) assuming Gaussian temper- 
ature fluctuations (Barreiro et al. 1997), to discriminate 
among standard CDM models. In this paper, we will use 
the correlation of excursion sets to test the Gaussianity of 
the CMB temperature fluctuations in the sky, including the 
presence of noise. An advantage of using this quantity is that 
it is very easy to identify all the excursion sets (i.e. pixels 
above a fixed threshold) in temperature maps. These correla- 
tions are amplified by a purely statistical effect with respect 
to those associated with the temperature field. They carry 
information on the 2-point probability density function, al- 
lowing us to distinguish between two different distributions 
even if the underlying 1-point probability density function 
is the same. The latter point will become clear with the use 
of some reference toy models. 

We discuss, in §2, the different 2-point probability den- 
sity functions generated in two different ways. §3 is dedicated 
to calculating the correlation of excursion sets for different 
models in the case of an ideal experiment whereas in §4 we 



Figure 1. The Gaussian (solid), x\ (short-dashed), X30 (dotted), 
Xgo (long-dashed) and Laplace (dotted-short dashed) l-pdf s with 
zero mean and unit variance are shown. 



include noise. Finally, in §5 and §6, we give the main results 
and conclusions of the paper, respectively. 



2 TWO POINT DISTRIBUTION FUNCTIONS 

We are interested in several Gaussian and non-Gaussian dis- 
tributions that could represent the temperature field pro- 
duced in different cosmological scenarios. Inflationary mod- 
els will generate Gaussian random fields whereas other mod- 
els such as those derived from topological defects will in 
general give rise to non- Gaussian ones. We shall introduce 
the chi-squared probability density function as a simple dis- 
tribution that asymptotically contains the Gaussian one. It 
has already been used in different cosmological contexts, in 
relation to the large scale structure (Coles 1989), to study 
topological properties of the CMB (Coles & Barrow 1987) 
and as a test of Gaussianity in the C0i?i5-DMR four-year 
sky maps (Kogut et al. 1995). 

The aim of this paper is to calculate the 2-point corre- 
lation of the excursion sets above a given threshold. In order 
to do this, we need to know the corresponding 2-point proba- 
bility density functions (hereinafter, 2-pdf's). The Gaussian, 
X^ with n degrees of freedom and Laplace 2-pdf's are ob- 
tained from the 1-pdf's (see figure 1) in two different ways: 
one will follow the "standard" procedure whereas the other 
represents a simple "toy" ansatz (Berry 1973, Jones 1996). 



2.1 Standard procedure 

a) Gaussian 

The 2-pdf for a Gaussian field with zero mean, variance 
a^ and correlation r is given by: 



fx{xi,X2) 



27ro-Vl--r2 



(1) 



-00 to 00. 



This distribution is symmetric and extends from 
b) Chi-squared Xn 

It is possible to generalize the univariate Xn distribution 
with n degrees of freedom to the bivariate case as follows: 
let us define the random vector y = {y 1,1/2) = (j'i{x\,x%) 
where Ox is a constant and x = (xi, 0:2) is a random variable 
described by a bivariate Gaussian satisfying (xi) = (a^i) = 
1, {x\X2) ~ Tx- Then, it is possible to calculate the 2-pdf for 
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y, as well as for the sum of n independent random vectors of 
this type. The special case r = then corresponds to a pair 
of independent Xn univariate distributions with n degrees of 
freedom multiplied by a^. In the Appendix, we derive the 
2-pdf for such a bivariate Xn distribution, resulting in the 
following expressions (cr^ = 2na^,T = r^): 
r < : 



fy2{yi,y2) = 



-yiy2n 



2cr2(l-r) V 2ro-' 



rn yi+y2 

'^ a(l-T) J„ 



1 



r(f) 

^/~2nTyiy2 



ail 



42(2/1,2/2) = 


1 /nyiy2\^ 
2/12/2 V 2a2 ) 
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[r(f)J 


T > : 




42(2/1,2/2) 


n f 1 


ny2n\- 

2ra^ J 


2cr2(l-r) I 




/TT B1+B2 
g V 2 ct(I-t) / 


^,(1 



(2) 



(3) 



r(f) 



(4) 



where Ji, and /^ are the Bessel and modified Bessel function 
of the first kind, respectively. 

The Xn distribution is asymmetric and extends from 
to 00. Since the temperature field is assumed to have zero 
mean, we will center this distribution in the following calcu- 
lations, which are also extended to negative values. In the 
limit n ^ cx), we recover a Gaussian process with mean /i, 



c) Laplace 

The Laplace 1-pdf is given by: 



fziz) = 



1 



-V2\i 



V2a 



(5) 



Following McGraw & Wagner (1968), it is possible to 
generalize the Laplace distribution to a 2-pdf if we assume el- 
liptical symmetry properties. An elliptically symmetric dis- 
tribution is characterized by a 2-pdf whose contour lines of 
equal height are ellipses in the (2:1,22) plane and satisfies: 



f{zi,Z2,T)=f(R,T) 



(6) 



where R — y 2^ -\- z^ — 2ziZ2r. The previous authors found 



for the Laplace 2-pdf: 
1 



/z(2l,22) = 



7ro"2\/l ~ T 



zA-Q 



y/2 {zf + 2| - 221 22r) 

crVl-r2 



(7) 



with mean value p = 0, variance a and correlation r. Ko is 
the modified Bessel function of second kind. This distribu- 
tion is symmetric and extends from — cxj to 00. Notice that 
for the case r = the Laplace 2-pdf is not given as the 
product of the 1-pdf's. 



2.2 Toy model procedure 

We will show in the present section a simple method to gen- 
eralize a 1-pdf to a 2-pdf (Berry 1973, Jones 1996). This 
generalization has been used by the first author to describe 



the statistical properties of echoes diffracted from rough sur- 
faces. If we have a process with correlation r and described 
by a 1-pdf f(x), with mean fi and variance a^, then, the 
distribution: 



fixi,X2) = f{xi)5{xi - X2)r + f{xi)f{x2){l - t) 



(8) 



is a 2-pdf with correlation r and marginal distributions given 
by the original 1-pdf. Here S(x — xo) is the Dirac delta dis- 
tribution centered on xq. This distribution extends in the 
same range as the f(x) being considered and satisfies the in- 
teresting properties: f{xi,X2) — f{xi)S{xi — X2) for r ^ 1 
and /(a;i,a;2) = f{xi)f{x2) for r — > 0, i.e., the same lim- 
its that can be obtained for the bivariate Gaussian. This 
generalization is a simple example which shows that knowl- 
edge of the 1-pdf and the correlation allows many possi- 
bilities for the 2-pdf. We will construct this 2-pdf for the 
Gaussian and non-Gaussian fields considered previously, and 
we will compare the correlation of regions obtained in this 
way with those obtained from the standard distributions 
given in §2.1. A CMB map associated to this toy model 
represents a surface consisting of fiat regions separated by 
vertical walls whose rms height is 2a, a being the disper- 
sion of the underlying 1-pdf (Berry, 1973). In particular, for 
the Gaussian 1-pdf this could mimic the patchy behaviour 
given by the Kaiser-Stebbins effect produced by a network 
of long straight strings (Kaiser & Stebbins 1984, Pompilio 
et al. 1995). For a simulation of the Kaiser-Stebbins effect 
see Magueijo & Lewin (1997). 



3 CORRELATION OF EXCURSION SETS 

Kaiser (1984) has calculated the 2-point correlation function 
of excursion sets above a threshold ua for a 3D Gaussian 
field of variance a^ and mean value zero. We will follow 
the same procedure for the distributions considered in this 
paper. 

The probability that a randomly chosen point lies above 
a certain level va is: 



Pi = 



f{x)dx 



(9) 



If we choose another point at distance r from the first, the 
probability that the field at both points takes a value ex- 
ceeding that threshold is: 



00 /"OC 

P2 = I / f{xi,X2,T)dxidX2 



(10) 



where t is the correlation of the field. 

Therefore, the 2-point correlation function for the ex- 
cursion sets, ^>y, is: 



i + C>40 = ^ 



P2_ 



(11) 



^>,/(r) gives the fractional excess probability that a point 
X2 lies above va, given that x\ also exceeds that threshold 
and \ x\ ~ X2\= r . 

Next we shall introduce the normalized correlation 
C>i/(r) associated to the characteristic function of an ex- 
cursion set h(x) defined by: 



© 0000 RAS, MNRAS 000, 000-000 



R.B.Barreiro, J.L. Sanz, E. Martinez-Gonzalez and J. Silk 



h{x) = 1 if f(x) > f 
h{x) — otherwise. 



(12) 



Then the correlation Cy,^{r) is the standard cor- 
relation coefficient of h(x): C'yv{r) = [{h{xi)h{x2)} — 
{h{xi)} {h{x2))]/[ahixi)ahix2)] with al{x) = (h'^ix)) - 
{h{x)) . It is straightforward to obtain the relation: 



C>Ar) = 



A 



1-Pi 



c>4^-; 



(13) 



We note that P2(r = 1) — Pi for the standard Gaus- 
sian, standard chi-squared and toy model distributions. We 
have also proved that this result is valid for high thresh- 
olds in the standard Laplace case, and it applies for the i/'s 
considered in the present paper. 

^>^(r) at zero lag (r — 1) only contains information 
about the 1-pdf, however C>^(0) — 1 for any distribution. 
Since the information from the 2-pdf is only encoded in the 
shape of the correlation of excursion sets and not in its nor- 
malization, it is better to use C>i/(r) to test the discrimi- 
nating power of the 2-pdf. 

The standard and toy model procedures will be consid- 
ered separately. 

3.1 Standard procedure 

a) Gaussian 

For a Gaussian field Pi is given by; 



P. 4.*^ 



(14) 



e ' dt is the complementary error 



where erfc(x) = -^ J" 
function. 

One of the integrals which appears in P2 can also be 
evaluated, resulting: 



i + C>. 



iv^: 



,1/2 



erfc^ 



dye 2 erfc 




in agreement with Kaiser (1984). This last expression must 
be evaluated numerically when estimating C>^. 
b) Non-Gaussian 

As explained in §2.1, the chi-squared distribution has 
mean value /i, different from zero, and it must be centered if 
we want this distribution to describe the temperature field. 
To take this into account, the integrals appearing in Pi and 
P2 have been evaluated between ji + va and 00. This pro- 
cedure is completely equivalent to centering the Xn distri- 
bution and taking the limits as shown in equations (9) and 
(10). Then, Pi is given by: 



Pi = l 



r(i) 



(16) 



where 7(0, x) — J^ e^^t^^^dt is the incomplete gamma func- 
tion. 

For the Laplace case, we have: 



Pi = l-^e"^" 



u <0 



Figure 2. We plot the correlation of regions above thresholds 
Us = 2.5 (upper figures) and Us = 3.5 (lower figures) versus the 
correlation of the cosmological signal, Ts, for f2 = 1, FWHM=10', 
three levels of noise (from left to right, (Tjv(lO') = (0, 1, 3) X 10~^) 
and different standard distributions: Gaussian(solid line), Xi 
(short-dashed). Coles approximation for Xi (dotted- long dashed, 
only in the plots with no noise), X30 (dotted), x^o (long-dashed) 
and Laplace (dotted-short dashed). 



Pi 



-V2v 



u> 



For both the Xn a^nd Laplace distributions P2, 
therefore Cy^, must be numerically evaluated. 



(17) 



and 



3.2 Toy model procedure 

For a 2-pdf generated as explained in §2.2, P2 — rPi+Pf (1 — 
r) and so the correlation of regions has a very simple form: 



C>.(r)=r(r) , 

where r is the correlation of the process. 



(18) 



In the top left plot of figure 2, the normalized correla- 
tions of excursion sets above f = 2.5 versus r for a Gaussian, 
Xn with n=l,30,60 degrees of freedom and Laplace standard 
distributions are shown. The different distributions go to 
zero when r goes to zero, except for the Laplace case due 
to the fact that its 2-pdf does not tend to the product of 
the 1-pdf's. We see that there are large quantitative differ- 
ences between the Gaussian and non-Gaussian distributions. 
As expected, when n increases the Xn correlation approaches 
the Gaussian one. For comparison, the C>i, for Xi calculated 
by Coles (1989) using the approximation of Politzer & Wise 
(1984) is also plotted. It is shown that this approximation 
only applies acceptably well for small r's, as was already 
pointed out by Coles. For high values of r the approxima- 
tion clearly overestimates the actual normalized correlation 
of regions for the Xi distribution. In the bottom left corner 
of figure 2, the same distributions are plotted for a thresh- 
old f = 3.5. We see that increasing the threshold amplifies 
the difference between the Gaussian and non-Gaussian dis- 
tributions. Notice that in the case of the toy model, Cyi, is 
given by the simple relation of equation (18). We emphasize 
that even if we have two fields with the same correlation and 
1-pdf (as the case for the standard and toy model procedure 
with the same marginal pdf), we can discriminate between 
them using the correlation of regions. 
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4 CORRELATION OF EXCURSION SETS 
INCLUDING NOISE 

In a real experiment, the presence of instrumental noise 
modifies the correlations calculated in §3. In addition, the 
result will depend on the amplitude of the signal and then 
we need to choose specific cosmological models. Whereas 
for the standard inflationary scenario with Gaussian tem- 
perature fluctuations, the radiation power spectrum can be 
determined very accurately once the matter content is spec- 
ified, this is not the case for other alternative scenarios such 
as topological defects. We will assume the same radiation 
power spectrum for all the Gaussian and non-Gaussian dis- 
tributions that we consider. We expect that the differences 
in a more realistic case, where the appropriate power spectra 
are known, are at least as large as the ones found under this 
assumption. In particular, we consider in the present section 
a fiat CDM model (baryon content ilf, = 0.05, Hubble con- 
stant h=0.5, cosmological constant A = 0) with adiabatic 
fluctuations and a Harrison-Zel'dovich primordial spectrum, 
kindly provided by N.Sugiyama. The Ce's have been nor- 
malized to the COBE 2-year maps (Cayon et al. f996; this 
normalization does not appreciably change with the 4-year 
data) being the signal dispersion as = 3.5 x fO~® for an an- 
tenna of FWHM=fO' and a smoothing of the same width. 
We also include the effect of noise that is assumed to be 
Gaussian white noise. Following standard observational pro- 
cedures, we have filtered signal plus noise with a Gaussian 
with the same width as the antenna. 

The angular correlation function for the cosmological 
signal with a Gaussian beam profile is given by: 



C{a,af) 



An ^-^ 

£ 



{2i + l)Cie 



-2l(l + l)a 



f Pe f cos a) 



and its normalized correlation TsI 

_ C{a,Of) 

at = C(0,o-/) , 



(19) 

(20) 
(21) 



where a is the Gaussian dispersion of the antenna (a/ = 
0.425 xFWHM). 

The correlation tm and the dispersion on for the filtered 
noise are very well approximated by (a/ ;$ 0.1 rad): 



TAT 



2 



At 



Ana) 



(22) 
(23) 



We choose different levels of noise, fixing the noise 
amplitude An = (0, 1.9, 17) x lO^'^^ in order to obtain 
crjv = (0,1,3) X lO"'"' after smoothing with a 10' FWHM 
Gaussian window. These levels of noise cover the range 
of sensitivities expected for future experiments (e.g. MAP, 
Planck Surveyor Mission). 

The presence of noise produces distortions in the distri- 
butions considered. The new distribution is given by a sum 
of a Gaussian process (coming from the noise) with disper- 
sion On plus the process that characterizes the signal with 
dispersion Os- The new 1-pdf and 2-pdf are given by: 



/(^) 



fN{z-y)fs{y)dy 



(24) 



f{zi: 



Z2 



fNizi - yi,Z2 - y2)fsiyi,y2)dyi,dy2 , (25) 



where /jv and fs refer to the distributions of the noise and 
the signal respectively. The integrals are evaluated over the 
range covered by the given signal distribution. The variance 
of the new process is given by cr^ = ctJ -I- a% and its corre- 
lation, for angles greater than the coherence angle ^^of the 
noise (6*^ ~ \f2af), can be approximated by the expression 
Tt — Tb/{1 + SNR~'^) where SNR = OsJon is the signal-to- 
noise ratio. We note that we are interested in temperature 
thresholds Vs in units of a^, which is related to the threshold 
vt, in units of at, by Vt = Vsl\J\ + SNR^'^. 

We will again consider both procedures separately. 



4.1 Standard procedure 

a) Gaussian 

For a signal described by a Gaussian distribution, the 
new process remains Gaussian with dispersion at and cor- 
relation Tt defined as before. Therefore, the functional form 
of the correlation of excursion sets is the same as the one 
without noise given by equations (14) and (15), where the 
threshold and correlation are now vt and rt, respectively. 

b) Non-Gaussian 

Pi and P2 are again given by equations (9) and (10). 
For Pi we have: 



■Pi = 2 / /«(j/)crfc[3(j/)]dj/ 



(26) 

l^-y) SNR)/V2 with 
and Laplace cases, re- 
spectively. For the Laplace distribution Pi can be evaluated 
analytically. 

On the other hand, if we assume tat ^ 1, some of the 
integrals which appear in P2 can be evaluated, resulting in : 



where g{y) = {vtVl + SNR? + (■ 
/i = \/^as and ^ — Q for the Xn 



/s(yi,J/2)erfc [g{yi)\ orfc [5(^2)] dyidy2 



(27) 



This approximation is valid for a ^ 6^ . To obtain a 
better estimation of the normalized correlation of regions 
for small angles when noise is present, we have evaluated P2 
as a power series considering terms up to second order in tn 
and then interpolated using its known value at r = 1 (again 
P2(r = l)=Pi). 

4.2 Toy model procedure 

In this case Pi and P2 are also given by equations (26) and 
(27) where fs(y) is the 1-pdf corresponding to the signal 
(Gaussian, Xn or Laplace) and fs{yi,y2) is the 2-pdf con- 
structed using the toy model procedure. In particular, as- 
suming rjv ^ 1 we obtain for P2: 



Ts 

T 



f4y)ieric[g{y)]fdy+{l 



rs)P? 



(28) 



In figure 2, Cyv versus Ts are plotted for different levels 
of noise and thresholds Us (for ^ = 1 and FWHM=10'), for 
the Gaussian, Xn with n=l,30,60 and Laplace distributions 
following the standard procedure. The main effect of the 
noise is to produce a rapid fall of the normalized correlation 
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Figure 3. In the four upper figures we plot the norinaUzcd cor- 
relation of regions for the toy model (thicker line, labelled t) 
and standard (thinner line, labelled s) 2-pdf's derived from the 
same 1-pdf. In the lower figures comparisons between correlation 
of regions for Gaussian and non-Gaussian distributions derived 
from standard 2-pdf's are shown. In all the plots, we have taken 
Us = 2.5, S7 = 1, FWHM=10' and no noise. The different lines 
correspond, in all the cases, to the correlation of regions derived 
from Gaussian (solid), xf (short-dashed), X30 (dotted), Xgo (loi^S" 
dashed) and Laplace (dotted-short dashed) distribution. 



of regions at small scales. However, it is still possible to see 
clear differences between Gaussian and non-Gaussian distri- 
butions, up to levels of SNR~1. As before, increasing the 
threshold Ps produces an amplification of the differences. 



5 RESULTS 

We have applied the calculations of the previous sections to 
open and fiat CDM models {Q. = 1,0.3,0.1, with 0-^(10') = 
3.5 X 10~*) and considered three different levels of noise 
(aiv(10') = (0,l,3) X 10-5). 

Besides the instrumental noise, our calculations are af- 
fected by statistical errors with whole sky coverage. Since 
we are considering very small angular scales, we expect that 
the cosmic variance will introduce a small uncertainty in our 
results. 

As we have already pointed out, increasing the thresh- 
old Vs , enhances the difference between the considered corre- 
lations. However, this also produces a rapid fall in the num- 
ber of excursion sets above Vs and, thus, higher statistical er- 
rors. On the other hand, increasing the width of the antenna 
improves the SNR but drastically decreases the number of 
spots above Va. This shows that we must find a compro- 
mise for choosing the optimal parameters for our analysis. 
In particular, we have taken an antenna FWHM=10' and a 
threshold Us — 2.5 to present our results. 

The normalized correlation of excursion sets above a 
threshold gives us different information from that of the 2- 
point correlation of the temperature fiuctuations, since we 
obtain additional information about the 2-pdf. For this rea- 
son, our method is of especial interest for discriminating 
between different distributions of the temperature fluctua- 
tions for a given power spectrum. In fact, as we have already 
pointed out, even if we have two fields described by the same 
1-pdf and correlation, we can show that it is possible to dis- 
criminate between them through the correlation of regions. 

In the upper part of figure 3, we plot pairs of standard 



Figure 4. We show the correlation of regions for flat (solid line) 
and open models with f2 = 0.1 (dashed line, in upper figures) and 
f! = 0.3 (dotted line, in lower figures) for different standard 2- 
pdf's (from left to right, Gaussian, x|q and Laplace). For all the 
figures, we have considered i/g = 2.5, FWHM=10' and no noise. 



and toy model distributions derived from the same 1-pdf for 
n = 1, FWHM^IO', Ua = 2.5 and no noise. The drawn er- 
rors have been estimated for all the distributions from high 
resolution Gaussian simulations of the sky. We expect them 
to be a good approximation to the real errors for the xio £^nd 
xio distributions. In the rest of the cases, even if the errors 
were underestimated, we don't expect our results to be very 
much affected since the difference between distributions is, 
generally, very clear. In addition, performing a simple Pois- 
sonian calculation of errors suggests that they are larger in 
the Gaussian case. Using a x^ test, we obtain a null result for 
the hypothesis that each pair of curves is derived from the 
same population, with a very high confidence level (;<; 99%). 
Considering open models, fl = 0.3,0.1, or introducing noise 
up to ajv(lO') = 3 X lO^'' does not significantly affect these 
results. 

On the other hand, we can compare the Gaussian 
and non-Gaussian standard distributions. For Q = 1, 
FWHM=10', i^a ~ 2.5 and no noise (see lower plots of fig- 
ure 3), we can distinguish between the Gaussian and non- 
Gaussian {Xn, n=l,30,60 and Laplace) standard distribu- 
tions with a confidence level <: 99%. Increasing the level 
of noise up to (jjv(lO') = 3 x 10-"" or changing the power 
spectrum does not appreciably modify these results. 

We can also try to discriminate between flat and open 
models for a given standard distribution through C>^. For 
i^s = 2.5, FWHM=10' and no noise, we can differentiate 
between the flat and open (fi = 0.1 or il = 0.3) models with 
a confldence level ~ 99%. In Figure 4, we have represented 
the Gaussian, xio ^^d Laplace distributions for the three 
considered values of Q, Vs ~ 2.5, FWHM=10' and no noise. 

Considering a higher threshold amplify the difference 
among the considered distributions. However it also pro- 
duces a rapid fall in the number of excursion sets above u 
and thus higher statistical errors, what leads to smaller con- 
fidence levels in the separation of the models. On the other 
hand, the smaller the threshold the closer C>^ for the dif- 
ferent distributions, since in the limit v — > —00 all the pixels 
are above u and the correlation of regions is zero for all the 
distributions. 
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6 CONCLUSIONS 

We have presented a method for testing the Gaussianity of 
the CMB temperature fluctuations in the sky based on the 
normahzed correlation function of excursion sets above a 
given threshold. It can be used to distinguish between stan- 
dard inflationary scenarios of generation of density fluctu- 
ations and those based on topological defects as well as to 
search for any trace of non-Gaussianity due to systematic 
errors or foregrounds. Such a technique was first introduced 
in cosmology by Kaiser (1984) to study the correlation of 
galaxy clusters in the context of biased scenarios of galaxy 
formation. 

As an application, we have constructed 2-pdf's from 
the Gaussian, Xn and Laplace 1-pdf's in two different ways. 
From these, we have obtained the normalized correlation of 
excursion sets including different levels of noise. This corre- 
lation contains additional information to that of the simple 
radiation power spectrum. 

Our main conclusion is that, using the correlation of 
excursion sets above high thresholds (e.g. Us — 2, 3) on sub- 
degree scales, it is possible to discriminate between differ- 
ent distributions even if the 1-pdf and correlation are the 
same. In particular, we are able to clearly distinguish be- 
tween the Gaussian and xio cases even in the presence of 
certain levels of noise within the range of sensitivities ex- 
pected for future experiments. Increasing the threshold am- 
plify the differences among the considered distributions but 
also the statistical errors, leading to smaller confidence lev- 
els. 

Finally, the correlation of regions can also be used to 
compare different il models though it is less efficient than 
directly using the radiation power spectrum. 
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APPENDIX A: THE 2-PDF FOR THE x% 
DISTRIBUTION 

We can find the 2-pdf for Xn in the following way. First, 
we construct the 2-pdf for Xi ~ ^^ j where X is a bivariate 
Gaussian process with zero mean. If fx{xi,X2) denotes its 
2-pdf then: 

fxliv^'V^) ^\ J \ fx {xi ^xi{yi),X2 = x2{y2)) , (Al) 

being | J \ the Jacobian of the transformation (a;i,a;2) -^ 
(yi,2/2), with 2/1 = xi, 2/2 = xi: 



fx^{yi,y2) = 



1 



yi+y2 



rx^yiy2 



cosh 



ol{l 



(A2) 



where a^. and Tx are the variance and correlation of the 
Gaussian process, respectively. 

Taking into account that the characteristic function of 
a given distribution is defined as the Fourier transform: 



0(ti,t2)= / / d2/id2/2/(yi,2/2)e'<'^^^+*^''^) 

— oo — oo 

we get for Xa'- 

1 



0^2(ti,i2) 



(A3) 



(A4) 



^1 - 2icrl{ti + t2) - 4(j|(l - T'i)Ut2 

On the other hand, the characteristic function of a sum 
of n independent processes is the product of their character- 
istic functions, so we have for Xn • 



<f>^2{tl,t2) 



(I>^2{ti,t2) 



(A5) 



This distribution has variance a^ = 2ncr^ , correlation t = t^ 
and nrean /i = na^ = \/^^- Then, we can write 2 in 
ternrs of a and r: 



<^y2 [HM)-- 



1 - « VT'^(*1 + *2) - l'^'(l - ^)*1*2 



(A6) 



Inverting eq.(A3) we obtain the 2-pdf for a chi-squared 
process with n degrees of freedom given by equations (2-4) 
of 52.1. 
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